APPENDIX X 



Program for Baselining, Normalizing Interpolating Then 

Calculati ng Spectral Overlap Integrals 

C This program has a non-standard DO WHILE loop 

INTEGER NPTS, NMAX, ROWS, ITER 
INTEGER EOF1, FLERR1, FLERR2 
INTEGER EOF2, FLERR3, FLERR4 
INTEGER EOF3, FLERR5, FLERR6 
INTEGER FLERR7 
INTEGER i j 

CHARACTER*30 fhamel, fhame2, fhame3, fhame4 
CHARACTER* 3 0 fhame5, fhame6, fhame7 
P ARAMETER(NMAX=3 500, LAMDA=601) 
REAL x,xxl (NMAX) ,yyl (NMAX),INTERV1 
REAL xx2 (NMAX) ,yy2 (NMAX), INTERV2 

REAL xx3 (NMAX), yy3 (NMAX),ENTERV3 

REAL yil (NMAX), yi2 (NMAX), yi3 (NMAX), yc (NMAX) 

REAL area 

CHARACTER* 1 SUBSTR, INITAR, LIGHT, INTMED 

FLERR1=0 

FLERR2=0 

FLERR3=0 

FLERR4-0 

FLERR5=0 

FLERR6=0 

FLERR7=0 

EOF 1=0 

EOF2=0 

EOF3=0 

INTERV1=0 

INTERV2=0 

INTERV3=0 

area=0 

write(* *) 'Do you wish to output intermediate files? (Y/N)' 
read(*,'(A)') INTMED 

write(*,*) 'Do you wish to process a substrate file? (Y/N)' 
read(*, '(A)') SUBSTR 

IF ((SUBSTR.EQ.'Y') .OR. (SUBSTR.EQ.'y')) THEN 
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ITER=0 

do 5 ITER=1, NMAX 

xx 1 (ITER)=0 

yyl (ITER)=0 

yil (ITER)=0 
5 continue 

write(*,*) 'Enter the name of the input substrate file:' 
read(V (A)') fhamel 

open 

(UNIT= 1 1 ,FrLE=fhamel,STATUS='OLDMOSTAT=FLERR 1 ,E 
RR=101) 

ROWS=0 

do while ((EOF1.EQ.0) .AND. (ROWS.LT.NMAX)) 
ROWS=ROWS+l 

Read (1 l,*,IOSTAT=EOFl) xxl (ROWS), yyl (ROWS) 

end do 

close (UNIT=11) 
NPTS=0 

IF(EOFl.NE.O) THEN 
NPTS-ROWS-1 

write(V(I4,lX,A12)') NPTS, 'points read/ 

ELSE 

NPTS-ROWS 

write(V(A28, 14, A12)') 'Too many data points! First', 
+ NMAX, ' points read . . . ' 

END IF 

call baseln(yyl,NPTS) 

IF ((INTMED.EQ. ' Y') .OR. (INTMED.EQ. 'y ')) THEN 

write(*,*) 'Enter the name of the output substrate file:' 

read(*,' (A)') fhame2 

open 

(TJNIT=12,FILE=mame2,STATUS=='NEW , ,IOSTAT=FLERR2, 
ERR=102) 

write(*,*) 'Writing data. . . ' 
END IF 
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do 12 L AMD A 

x= (i-l)+200 

call locate(xxl,NPTS,xj) 

if ((j.eq.O). OR. (jeqNPTS)) then 
INTER V 1=0 

else 

intervi= ((yylO+l)-yyl(j)) / (xxifl+i) - xxiQ))) 
(x-xxl(j)) 

+ +yylG) 
end if 

yil (i)=INTERVl 

IF ((INTMED.EQ. ' Y') .OR. (INTMED.EQ. 'y ')) THEN 

if (O'.eq.O) .OR. G'.eq.NPTS)) then 

GO TO 12 

else 

write(12,'(lx, f7.2,i6 > 3fl2.2)')xj,xxl(j),xxl0+l), 
+ INTERV1 

endif 

END IF 

12 continue 

IF ((INTMED.EQ. 'Y') .OR. (INTMED.EQ. 'y')) THEN 
close (UNIT=12) 

END IF 

ELSE 

ITER=0 

do 14 ITER^l, LAMDA 
yil (ITER) =1 
14 continue 

ENDIF 
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write(*,*) 'Do you wish to process an initiator file? (Y/N)' 
read(V(A)') INITAR 

IF ((INITAR.EQ.'Y') .OR. (INITAR.EQ. V)) THEN 
ITER=0 

do 15ITER=1,NMAX 

xx2 (ITER) =0 

yy2 (ITER)=0 

yi2 (ITER) =0 
1 5 continue 

write(*,*) 'Enter the name of the initiator file:' 
read(V(A)') mame3 

open 

(UNIT= 1 3 ,FILE=mame3 ,STATUS='OLD\IOSTAT=FLERR3, 
ERR- 103) 

ROWS=0 

do while ((EOF2.EQ.0) .AND. (ROWS.LT.NMAX)) 
ROWS=ROWS+l 

read(13,*,IOSTAT=EOF2) xx2 (ROWS), yy2 (ROWS) 

end do 

close (UNIT- 13) 
NPTS=0 

IF (EOF2.NE.0) THEN 
NPTS=ROWS-l 

write(*,l(14,lX,A12)') NPTS,'points read.' 

ELSE 

NPTS=ROWS 

write (V(A28,I4,A12)') 'Too many data points! First \ 
+ NMAX,' points read...' 

END IF 

call baseln (yy2,NPTS) 

IF ((INTMED.EQ.'Y').OR. (INTMED.EQ.'y')) THEN 

write(*,*) 'Enter the name of the output initiator file: ' 
read( V(A)') fhame4 
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open 

(UNIT= 1 4,FILE=fhame4,STATUS= 'NEW ' ,IOST AT=FLERR4, 
ERR=104) 

vvrite(*,*) 4 Writing data. . . ' 

END IF 

x=0 

i=0 

do 22 i=l,LAMDA 
x= (i-l)+200 

call locate (xx2,NPTS,x,j) 
if ((j.eq.O).OR. (j.eq.NPTS)) then 
INTERV2-0 

else 

INTERV2= ((yy2(j+l)-yy2(j)) / (xx2(j+l) - xx2 
0))) * (x-xx2 (j)) 
+ +yy2 (j) 

end if 

yi2 (i)=INTERV2 

IF ((INTMED.EQ. *Y*) .OR. (INTMED.EQ-. y)) THEN 

if ((j.eq.O) .OR. (j.eq.NPTS)) then 
GO TO 22 

else 

write(14,'(lx,f7.2,i6,3fl2.2)') x, j, xx2 

G),xx2 0+1), 

+ INTERV2 
endif 
END IF 
22 continue 

IF ((INTMED.EQ. 'Y').OR. (INTMED.EQ. y )) THEN 

close (UNIT- 14) 

END IF 
ELSE 

ITER=0 
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do 24 ITER=1,LAMDA 

yi2 (ITER)=1 
24 continue 



ENDIF 

write(*,*) 'Do you wish to process a light source file? (Y/N)' 
read(*,'(A)') LIGHT 

IF ((LIGHT.EQ.'Y').OR. (LIGHT.EQ.'y')) THEN 
ITER=0 

do 25 ITER=1,NMAX 
xx3 (ITER) =0 
yy3 (ITER) =0 
yi3 (ITER)=0 
25 continue 

write(*,*) 'Enter the name of the light source file:' 
read(*,' (A)') fhame5 

open 

(UNIT=15,FILE=fhame5,STATUS='OLD',IOSTAT=FLERR5, 
ERR- 105) 

ROWS=0 

do while ((EOF3.EQ.0) .AND. (ROWS.LT.NMAX)) 
ROWS=ROWS+l 

read(15,*JOSTAT=EOF3) xx3 (ROWS), yy3 (ROWS) 

end do 

close (UNIT= 15) 
NPTS=0 

IF (EOF3.NE.0) THEN 
NPTS=ROWS-l 

write(*,l (14,1X,A12)1)NPTS, 'points read. ' 

ELSE 

NPTS=ROWS 

write (*,' (A28,I4,A12)') 'Too many data points! First 
+ NMAX, ' points read 

END IF 

call norm (yy3,NPTS) 
IF ((INTMED.EQ.'Y').OR. (INTMED.EQ. 'y ')) THEN 

write(*,*) 'Enter the name of the light source output file:' 
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read(V (A)») fhame6 
open 

(UNIT=16,FILE-fname6,STATUS- , NEW\IOSTAT=FLERR6, 
ERR=106) 

write(*,*) 'Writing data...' 

END IF 
x=0 
i=0 
j-0 

do 32 i=l,LAMDA 
x- (i-l)+200 

call locate (xx3,NPTS,xj) 
if (0-eq.O) .OR. (j.eq.NPTS)) then 
INTERV3-0 

else 

INTERV3= ((yy3(j+l) - yy3fl)) / (xx3(j+l) - 
xx3(j))) * (x-xx3 (j)) 

+ +yy3 0) 

end if 

yi3 (i) =INTERV3 

IF ((INTMED.EQ. ' Y') .OR. (INTMED.EQ. 'y')) THEN 

if (0'.eq.0).OR. O'.eq.NPTS)) then 
GO TO 32 

else 

write(16/(lx,f7.2,i6,3fl2.2)') xj,xx3 (j),xx3 0+1), 
+ INTER V3 

endif 

END IF 

32 continue 

IF ((INTMED.EQ/ Y') .OR. (INTMED.EQ. 'y')) THEN 

close (UNIT= 16) 

END IF 

ELSE 
ITER=0 

do 34 ITER=1,LAMDA 

yi3 (ITER) -1 
34 continue 
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END IF 



ITER=0 

DO 40 ITER= 1 ,L AMD A 

yc (ITER)=0 
40 CONTINUE 

DO 55 i=l,LAMDA 

yc (i)-yil (i)*yi2 (i)*yi3 (i) 
55 CONTINUE 

write(*,*) 'Enter the filename for cumulative data:' 

read (* '(A)') fhame7 

open 

0JNIT=17,FILE=mame7,STATUS= , NEW , ,IOSTAT=FLERR7, 
ERR=107) 

write(*,*) 'Writing data. . . ' 
CALL integ (yc,LAMDA,area) 

write(*,'(lX,A26,All,F12.6)') 'The area under the product', 
+' curve is:', area 

write(17,'(lX,A26,Al 1,F12.6)') 'The area under the product', 
+' curve is: ',area 

DO 60 i=l,LAMDA 
x= (i-l)+200 

write(17,' (1X,F6.1,2F11.2,F11.4,E15.6)') x,yil (i),yi2 

(0, 

+ yi3 (i) ,yc (i) 

60 CONTINUE 

close (UNIT- 17) 

101 IF (FLERR1 .NE. 0) THEN 
write(*,*) 'Unable to open substrate file!' 

END IF 

102 IF (FLERR2 .NE. 0) THEN 

write(*,*) 'Unable to create substrate output file!' 
END IF 

1 03 IF (FLERR3 .NE. 0) THEN 
write(*,*) 'Unable to open initiator file!' 



31 



END IF 



1 04 IF (FLERR4 .NE. 0) THEN 

write(*,*) 'Unable to create initiator output file!' 
END IF 

1 05 IF (FLERR6 .NE. 0) THEN 

write(*,*) 'Unable to open light source file!' 
END IF 

1 06 IF (FLERR6 .NE. 0) THEN 

write(*,*) 'Unable to create light source output file!* 
END IF 

107 IF (FLERR7 .NE. 0) THEN 

write(*,*) 'Unable to create cumulative output file!* 
END IF 

write(*,*) 'Program exiting normally. . . ' 
END 

SUBROUTINE locate (xx,n,xj) 

INTEGER j,n 

REAL x,xx (n) 

INTEGER j 1 jm ju 

J1=0 

ju=Tl+l 

10 if G'u-jl.gt.l) then 

jm=0u+jl)/2 

if ((xx (n).ge.xx (1)) .eqv. (x.ge.xx Om))) then 
jl=jm 

else 

ju-jm 

endif 
goto 10 
endif 

if (x.eq.xx (l))then 
j=l 

else if (x.eq.xx(n))then 
else 

j=jl 

endif 

return END 

SUBROUTINE baseln (yy,N) 
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INTEGER N, i 

REAL yy (N), minno, temp 

minno=yy ( 1 ) 
i=0 

temp=0 

DO 10i=2,N 

IF (yy (i) XT. minno) THEN 
minno=yy (i) 

END IF 
10 CONTINUE 

DO 20 i=l,N 

temp=yy (i)-minno 

yy (i)=temp 
20 CONTINUE 

END 

SUBROUTINE norm (yy,N) 

INTEGER N,i 

REAL yy (N),maxno,temp 

maxno=yy ( 1 ) 
i=0 

temp=0 

DO 5 i=2,N 

IF (yy(i).GT.maxno) THEN 
maxno=yy (i) 

END IF 
5 CONTINUE 

i=0 

DO 10, i=l,N 

temp=yy (i)/maxno 

yy (i)=temp 
10 CONTINUE 

END 
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SUBROUTINE integ (yy,N,area) 

INTEGER N,i 

REAL yy (N), sum, area 

i=0 

sum=0 
area=0 

DO 10 i=l,N-l 

sum=sum+ (yy(i)+yy (i+l))*0.5 
10 CONTINUE 

area=sum/ 100000 

END 
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APPENDIX 2 

Program to create uniformly spaced csv data from unevenly spaced 
tabular data 

#include <stdio.h> 
#include <stdlib.h> 
#include <math.h> 

#defme NMAX 3501 
#define STRMAX 151 
#define FNMAX 81 
#defme OUTPTS 801 

void Iocate(float xx[], unsigned long n, float x, unsigned long *j); 
void norm(float xx[], unsigned long int n); 
void baseline(float xx[], unsigned long int n); 

int main() 
{ 

char fiiamein(FNMAX], string[STRMAX], *str, ptr, 
fiiameout[FNMAX], another; 

float xdata(NMAX], ydata[NMAX], xinter[OUTPTS+l], 

yinter[OUTPTS+l]; 

unsigned long int index, i, j; 

int choice; 

FILE *fpin, *fpout; 

another = *Y'; 

do { 

for (i = 0; i <= NMAX-1; i++) { 
xdata [i] = 0; 

ydata [i] = 0; 

}; 

for (i - 0; i <= OUTPTS; i++) { 
xinter [i] = 0; 

yinter (i] = 0; 

}; 

printfC'Enter name of the input file (80 chars max, no spaces): M ); 
scanf("%s", fhamein); 

printf("File name is %s\n ,, ,friamein); 

fpin = fopen(fhamein, M r"); 

if(fpin = NULL) { 
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printf( M Cannot open %s\n",fhamein); 
exit(l); 

}; 

index = 1 ; 

while (1) { 

str_ptr = fgets(string,STRMAX-l,fpin); 
if(str_ptr = NULL) 

break; 
if(index = NMAX) 

break; 

sscanf (string, "%f %f ' , &xdata [index] , &ydata [index]); 
index++; 

}; 

f close (fpin); 

if((index = NMAX) && (str_ptr != NULL)) { 
index—; 

printf('Too many data points! Using first %d points 
only... \n", index); 
} 

else { 

index—; 

printf("%d points read... \n", index); 

}; 

printf("\nEnter option for data processing\n"); 
printf(" 1 : Normalize the data after interpolation^"); 
printf("2: Baseline the data after interpolation^"); 
printf("3: First interpolate, then baseline and finally "); 
printf("normalize the data\n"); 
printf( M 4: Simply interpolate the data\n M ); 
printf("5: Simply normalize the data\n"); 
printf("6: Simply baseline the data\n M ); 
printf("or\n"); 

printf("0: to exit the program without any data processing\n"); 

printf( M \nEnter option (0-6): "); 

scanf("%d M ,&choice); 

if (choice = 0) 
exit(2); 

printf("\nEnter name of the output file (80 chars max, no spaces): "); 
scanf( M %s",fhameout); 

printf( M File name is %s\n",mameout); 
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fpout = fopen(fhameout,"w"); 

if(fpout = NULL) { 

printf(" Cannot open %s\n",fhameout); 
exit(3); 

}; 

for(i = 1 ; i <= OUTPTS; i++) { 
xinter[i] = 200+((float)i-l); 
locate(xdata,index,xinter[i],&j); 
if(G = 0)||G = index)) 
yinter[i] = 0; 

else 

yinter [i] = (xinter [i] -xdata [j ]) * ((ydata 0+1] -ydata [j ]) / 
(xdata [j+1] -xdata [j ])) +ydata [j ] ; 

}; 

if ((choice = 2) |[ (choice = 3)) { 
baseline(yinter,OUTPTS); 

}; 

if (choice = 6) { 

baseline(ydata,index); 

}; 

if ((choice = 1) || (choice = 3)) { 
norm(yinter,OUTPTS); 

}; 

if (choice = 5) { 

norm(ydata, index); 

}; 

if ((choice >= 1) && (choice <= 4)) { 
for (i = 1 ; i <= OUTPTS- 1 ; i++) { 
fprintf(fpout,"%13.5E, ",yinter[i]); 

}; 

fprintf(fjpout, ,,0 /ol3.5E\n M ,yinter[OUTPTS]); 
} 

else 

if ((choice = 5) || (choice = 6)) { 
for (i = 1 ; i <= index- 1 ; i++) 
fprintf(fpout,"%13.5E, ",ydata[i]); 

}; 

fprintf(fpout,"%13.5E\n",ydata[index]); 

}; 

fclose( fpout); 
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printf("File %s written.\n\n",fhameout); 
printf("Process another file (Y/y/N/n)?: "); 
scanf("%ls",&another); 
} while (another = 'Y' || another = '/); 

printf("Exiting...\n"); 
return(O); 

} 

void locate(float xx[], unsigned long n, float x, unsigned long *j) 

{ 

unsigned long jujmjl; 
int ascnd; 

jl=0; 

ju=n+l; 

ascnd=(xx[n] >= xx[l]); 
while (ju-jl> 1) { 
jm=(ju+j I) » 1; 
if (x >= xx [j m] = ascnd) 

jl=jm; 
else 

ju=jm; 

} 

if(x = xx [1] ) 
else if (x = xx[n]) 

*j=n-i; 

else 

*j=ji; 

} 

void norm(float xx[], unsigned long int n) 

{ 

unsigned long int i; 
float maxdata, temp; 

maxdata = xx [1] ; 
temp - 0; 

for(i = 2; i <= n; i++) { 
if(xx[i] > maxdata) 
maxdata = xx[i]; 

}; 

for(i = 1 ; i <= n; i++) { 
temp = xx[i]/maxdata; 
xx [i] = temp; 

}; 



38 



} 



void baseline(float xx[], unsigned long int n) 
{ 

unsigned long int i; 
float mindata, temp; 
mindata = xx [ 1 ] ; 
temp = 0; 



10 for(i - 2; i <= n; i++) { 

if(xx[i] < mindata) 
mindata = xx [i] ; 

}; 

15 for(i = l;i<=n;i++) { 

temp = xx(i] - mindata; 
xx(i] = temp; 

}; 

M } 

Q 

□ 



APPENDIX 3 



Program for Determining Strength of Wavelength Response in a 

Region 

#include 

<stdio.h> 

#include 

<stdlib.h> 

#include 

<math.h> 

#define NMAX 
3501 

#define STRMAX 
151 

#define FNMAX 
81 

#define OUTPTS 
801 

void locate(float xxQ, unsigned long n, float x, 
unsigned long *j); 

void norm(float xx[], unsigned long 
int n); 

void baseline(float xxQ, unsigned 
long int n); 

void partinteg(fioat xxQ, unsigned long intxl, unsigned long int 

x2, 

float *area); 
int 

main() 
{ 

charfnameinfFNMAX], string[STRMAX], *str_ptr, 
fnameoutfFNMAX], another; 

float xdata[NMAX], ydata[NMAX], xinter[OUTPTS+1], 
yinter[OUTPTS+1]; 

float 

totalarea.aA.aB.aC.aD.aE.aF.aG.aH.al.aJ.a 

K; 

unsigned long int index, i, 

j; 

int choice; 
FILE *fpin, 
*fpout; 
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another = 'V; 

printff'Contact Rajdeep S. Kalgutkar, SRC-CRC 7-3003, for 
further info\n H ); 

do{ 

for (i = 0; i<= NMAX-1; 
i++){ 

xdata[i]=0; 
ydata[i]=0; 

}; 

for (i = 0; i <= OUTPTS; 
i ++ ){ 

xinter[i]=0; 
yinter[i]=0; 

}; 



printf("\n Enter name of the input file (80 chars max, no 
spaces): "); 

scanf("%s",fnamein); 



fpin = 
fopen(fnamein,"r"); 

if (fpin == 
NULL) { 

printffCannot open %s. 
Exiting. ..\n",fnamein); 

exit(1); 

}; 

index = 1 ; 

while (1) { 
str_ptr = fgets(string,STRMAX- 
1 ,fpin); 

if(str_ptr == NULL) 

break; 
if(index == NMAX) 

break; 
sscanf(string,"%f 
%f',&xdata[index],&ydata[index]); 



printf("File name is 
%s\n",fnamein); 
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index++; 

}; 

fclose(fpin); 

if((index == NMAX) && (str_ptr != NULL)) { 
index—; 

printf("Too many data points! Using first %d points 
only...\n",index); 
} 

else { 
index—; 

printf("%d points 
read...\n",index); 

}; 

printf("\n Enter option for data 
^ processing\n n ); 

printf("1 : Simply interpolate the 
p! data\n"); 

V printf("2: Normalize the data after 

uj interpolation^"); 

<0 printf("3: Baseline the data after 

0 interpolation\n"); 

1 printf("4: First interpolate, then baseline and finally 

printf("normalize the 
ry data\n M ); 
D printf("or\n"); 

H printf("0: to exit the program without any data 

processing\n"); 

printf("\nEnter option (0- 
4): "); 

scanf("%d",&choice); 
if (choice == 0) 
exit(2); 

printf("\nEnter name of the output file (80 chars max, no 
spaces): "); 

scanf("%s ,, I fnameout); 

printf("File name is 
%s\n",fnameout); 

fpout = 
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fopen(fnameout,"w"); 

if (fpout == 
NULL){ 

printf("Cannot open %s. 
Exiting.. An", fnameout); 

exit(3); 

}; 

for(i = 1 ; i <= OUTPTS; 
i++){ 

xinter[i] = 200+((float)i- 

D; 

locate(xdata,index,xinter[i],&j); 

if(0 == 0) || 0 == 
index)) 

yinterfi] = 0; 
else 

yinter[i]=(xinter[i]-xdata[j])*((ydata[j+1]- 
ydata[j])/ 

(xdataO+1 ]-xdata[j]))+ydata[j]; 

}; 

if ((choice == 3) || (choice == 4)) { 
baseline(yinter,OUTPTS); 

}; 

if ((choice == 2) || (choice == 4)) { 
norm(yinter.OUTPTS); 

}; 



partinteg(yinter,51,OUTPTS,&totala 
rea); 



partinteg(yinter,51,101,&a 
A); 

partinteg(yinter,101 ,1 51 ,& 
aB); 

partinteg(yinter,1 51 ,201 ,& 
aC); 
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partinteg(yinter,201 ,251 ,& 
aD); 

partinteg(yinter,251 ,301 ,& 
aE); 

partinteg(yinter,301 ,351 ,& 
aF); 

partinteg(yinter,351 ,401 ,& 
aG); 

partinteg(yinter,401 ,451 ,& 
aH); 

partinteg(yinter,451 ,501 ,& 



n partinteg(yinter,501 ,551 ,& 

t aJ); 

UJ partinteg(yinter,551 ,OUTPTS,&aK); 

u fprintf(fpout,"The total area is: 

JL L %14.6E\n",totalarea); 
p fprintf(fpout,'The area under region A is: 

fU %6.2f%%\n",aA*1 00/totalarea); 

Rj fprintf(fpout,"The area under region B is: 

Q %6.2f%%\n",aB*1 00/totalarea); 

"^■ s fprintf(fpout,"The area under region C is: 

%6.2f%%\n",aC*1 00/totalarea); 

fprintf(fpout,"The area under region D is: 
%6.2f%%\n",aD*1 00/totalarea); 

fprintf(fpout, M The area under region E is: 
%6.2f%%\n",aE*1 00/totalarea); 

fprintf(fpout,'The area under region F is: 
%6.2f%%\n",aF*1 00/totalarea); 

fprintf(fpout,"The area under region G is: 
%6.2f%%\n",aG*1 00/totalarea); 

fprintf(fpout,"The area under region H is: 
%6.2f%%\n",aH*1 00/totalarea); 

fprintf(fpout,"The area under region I is: 
%6.2f%%\n",al*1 00/totalarea); 

fprintf(fpout,"The area under region J is: 
%6.2f%%\n",aJ*1 00/totalarea); 

fprintf(fpout, M The area under region K is: 
%6.2f%%\n\n H ,aK*100/totalarea); 
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for (i = 1 ; i <= OUTPTS-1 ; { 
fprintf(fpout,"%13.5E, ",yinter[i]); 

}; 

fprintf(fpout,"%13.5E\n ,, ,yinter[i]); 
fclose(fpout); 

printf("File %s 
written.\n\n*\fnameout); 

printf("Process another file 
(Y/y/N/n)?: "); 

scanf("%1s",&another); 
} while (another == "Y 1 1| another 
== V); 



printf( M Exiting..An 

"); 

return(0 

); 

} 

void locate(float xxQ, unsigned long n, float 

unsigned long *j) 

{ 

unsigned long jujmjl; 
int ascnd; 

jl=0; 

ju=n+1; 

ascnd=(xx[n] >= xx[1]); 
while (ju-jl > 1) { 

jm=(ju+jl) » 1; 

if (x >= xx[jm] == ascnd) 

jl=jm; 
else 

ju=jm; 
} 

if (x == xx[1]) 
1=1; 
else if(x == 
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xx[n]) 
*j=n- 

1; 

else 
*j=jl; 

} 

void norm(float xxQ, unsigned long 

int n) 

{ 

unsigned long 
int i; 

float maxdata, 
temp; 

maxdata = 
xx[1]; 
temp = 

0; 

for(i = 2; i <= n; { 
if(xx[i] > 
maxdata) 

maxdata = 
xx[i]; 

}; 

for(i = 1 ; i <= n; i++) { 
temp = xx[i]/maxdata; 
xx[i] = temp; 

}; 

} 

void baseline(float xx[|, unsigned 
long int n) 

{ 

unsigned long 
int i; 

float mindata, 
temp; 

mindata = xx[1]; 
temp = 

0; 

for(i = 2; i <= n; { 
if(xx[i] < 
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mindata) 

mindata = 
xx[i]; 

}; 

for(i = 1; i <= n; { 
temp = xx[i] - mindata; 
xx[i] = temp; 

}; 

} 

void partinteg(float xx[], unsigned long int x1, unsigned long int 
x2, 
float *area) 

{ 

unsigned long 
int i; 

float temp; 

temp = 

0; 

for(i = x1; i <=x2 - 1; 

{ 

temp = temp + (xx[i] + xx[i+1])/2; 

}; 

*area = temp; 

} 
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APPENDIX 4 



SRC Curing Resource dB 4 Query Select2 

Sub Initialize 

Dim ses ses As New NotesSession 
Dim db jdb As NotesDatabase 
Dim view view As Notes View 

Dim note _notel As NotesDocument, note_note2 As NotesDocument 
Dim i cnt As Integer, i add As Integer 

Set db db =ses ses.CurrentDatabase 

Set note notel=ses ses.DocumentContext 

Redim Preserve arr WavelengthRegion(O) air WavelengthRegion(O) _ 

If note_notel.Selectionl(0) o Or rwte_notel.Selection2(0) o Then 
If note notel.Selectionl(O) o Then 
Set view_view = db_db.GetView("By NotelD") 

If note notel.Selection2(0) o '"Then 

Set note_note2=view_view.GetDocumentByKey(Right("00000000" & 

note_notel.Selection2(01, 8)) 

Else 

Set note note2 = view_view. GetDocumentBy Key (Right( "00000000" & 
note_notel.Selectionl(0), $)) 
End If 

If Not (note_note21s Nothing) Then 

If note note2.Hasltem('WavelengthRegion") Then 

i cnt = - 1 

Forall vals In note_note2.WavelengthRegion 
If valso""Then 

i_cnt = i_cnt + 1 

Redim Preserve arr WavelengthRegion(i_cnt) 
arr_WavelengthRegion(i cnt) = vals 

End If 
End Forall 

End If 

End If End If 

If note notel.Type(O) = "S" Then 

Set view- View = db db.GetView("Substrate") 
Elseif note _notel.Type(0) = "P" Then 
Set view-view = db_db.GetView("InitiatorSensitizer") Else 
Set view-view = db_db.GetView("LightSource") End If 
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* Set note note2 =view-view. GetFirstDocument 
i_cnt = - 1 

Do While Not (note note2 Is Nothing) 
If note_note2.Name(0) o Then 
L add = True 

If arr WavelengthRegion(O) o Then 
Ladd = False 

Forall valsl In note _note2.WavelengthRegion 
Forall vals2 In arr_WavelengthRegion 
If valsl = vals2 Then 
iadd = True 
Exit Forall 

End If 
End Forall 

If Ladd Then 

Exit Forall 

End If 
End Forall 
End If 

If L add Then 

i cnt = i cnt + 1 

Redim Preserve arr_names(i cnt) 

arr names(i_cnt) = note_note2.Name(0) 

End If 

End It 

Set note note2 = view-view. GetNextDocument(note_note2) 

Loop 

note _notel .Names = arr_names 

End Sub 

SRC Curing Resource dB 4 Query Select2 Save Agent 

Sub Initialize 

Dim ses sesAs New NotesSession 
Dim db_db As NotesDatabase 
Dim view view As NotesView 

Dim note_notel As NotesDocument, note_note2 As NotesDocument 

Set db_db = ses ses.CurrentDatabase 

Set note note) = ses ses.DocumentContext 
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Select Case note_notel.Type(0) 
Case"S" 

Set view view = db_db. Get Vie w(' '(Substrate)") 

Set note_note2=view view.GetDocumentByKey(note_noteLSubstrate(0)) 
Case "P" 

Set view-view db-db.GetView("(InitiatorSensitizer)") 

Set note note2 -view view.GetDocumentByKey(note notel.Photolnitiator(O)) 
Case "L" 

Set view-view =db db.GetView("(LightSource)") 

Set note_note2 = view_view.GetDocumentByKey(note_noteLLightSource(0)) 
End Select 

If note_notel.MexWction(0) = "Add" Then 
If note_notel.Selectionl (0) o *'" Then 

Print "[!" + note notel.dbname(O) + "/Query Selectionl?OpenForm&" & 

note_notel.Selectionl (0) & "&" & note_note2.Noteld & ")" 

Else 

Print "[!" + note notel.dbname(O) + "/Query Selectionl?OpenForm&" & 
note_note2.Noteld & "]" 

End If Elseif note notel.NextAction(O) = "Separate" Then 
If note_notel.Selectionl(0) o "" Then 

Print "[/" + note notel.dbname(O) + "/QuerySelectionl?OpenForm&" & 

note_notel.Selectionl(0) & "&" & note note2.Noteld & ")" 

Else 

Print "[!" + note notel.dbname(O) + "/Query Selectionl?OpenForm&" & 
note_note2.Noteld & ")" 
End If 

Elseif note _notel.NextAction(0) = "Separate" Then 
If note _notel.Selection2(0) o "" Then 

Print "[/" + note notel.dbname(O) + "/QuerySelectionResults?OpenForm&" & 

note_notel. Selection) (0) & "&" & note note).Selection2(0) & 

"&" & note_note2.Noteld & "]" 

Elseif note notel.Selectionl(O) o "" Then 

Print- "[" + note notel.dbname(O) + < VQuerySelectionResults?OpenForm&" & 

note_notel.Selectionl (0) & "&" & note note2.Noteld & ")" 

Else 

Print "[/" + note_notel.dbname(0) + "lQuerySelectionResults?OpenForm&" & 
note_noteMoteld & "]" 
End If 
Else 

If note notel.Selection2(0) o "" Then 
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Print "[/" + note_notel.dbname(0) + "/QuerySelecnonOverlayResults?OpenForm&" & 
note_notel.Selectionl(0) & & note_notel.Selection2(0) & 

& note_note2.Noteld & ")" 
Elseif notenotel.Selectionl(O) o ""Then 

Print "[/" + note_notel.dbname(0) + "/QuerySelecnonOverlayResults?OpenForm&" 

& note_notel. Selection) (0) & "&" & note note2.Moteld & ")" 

Else 

Print "[/" + note_notel.dbname(0) + "lQuerySelectionOverlayResults?OpenForm&" 
& note_note2 .Noteld & ")" 
End If 

End If 
End Sub 
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SRC Curing Resource dB 4 Query Overlay Open Agent 

Sub Initialize 

pirn ses_ses As New NotesSession 
Dim db db As NotesDatabase 
Dim view-view As NotesView 

Dim note notel As NotesDocument. note__note2As NotesDocument 
Dim i_cntAs Integer, i addAs Integer 

Set db db = ses ses.CurrentDatabase 

Set note notel =ses ses.DocumentContext 

S et view- view = db_db. G etView("By N otel D ") 



If note notel. Selectionl(O) o "" Then 

S et note-note2 = view-view. G etD ocumentByKey(R fight("00000000" + notenotel. 
M S election"! (0). 8)) 

□ 

O 20 If Not (note _note2 Is Nothing) Then 

note _notel.datal = note_note2.EmissData 
note notel. maxfreql = note_note2.MaxFreq 
Jn End If End If 

5 

■ 25 If note_notel.Selection2(0) o Then 

H Set note-note2 = view view.GetDocumentByKey(Right("00000000" + 

note_notel.Selection2(0), 8)) 



^ If Not (note note2 Is Nothing) Then 

P0 note notel.data2=note_note2.EmissData 

note _notel.maxfreq2 = note_note2.MaxFreq 
End If End If 

If note notel. Selection3(0) o Then 
35 Set note-note2 = view_view.GetDocumentByKey(Right("00000000" + 
note_notel.Selection3(0), 8)) 

If Not (note _note21s Nothing) Then 

note notel. data3 = note note2.EmissData 
40 note _notel.maxfreq3=note note2.MaxFreq 

End If End If End Sub 
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* » 

APPENDIX 5 

import Java. awt.*; import java.awt. event.*; import java.applet.*; 
5 public class SRC Charts extends Applet { int gi count; 

double GetHMax(String str_in) { String str_current; double dbl hmax; 

str_current = ""; dblhmax = 0; for Tint i cnt = 0; i_cnt < str_in.length(); i cnt++) { 
10 if(str in. region Match es(i cnt, 0, 2)) 

if(Double.valueOf(str current). double Value() > dbl Jimax) 

dbl hmax = Double. valueOf(str current).doubleValue(); 

str current - *" ; 

i_cnt++; 
1 5 gi count++; ) else { 

str current = str current. concat(str in. substring^ cnt, i cnt + 1)); 

g ) 

g if(str_current.length() > 0) { 

S 20 if(Double. valueOf(str_current).doubleValue() > dbl_hmax) 

jg dbl hmax « Double. valueOf(str_current).doubleValue(); gi count-H-; ) return dbl hmax; 

□ int StringTolnt( String str_in, double dbl hmax) { double dbl_pos; 

l_25 

p dbl_pos - getSize().height - (25 + (Double. value0f(str_in). double ValueO * 

m ((getSizeO-height - 50) / dbl hmax))); return (int)dbl-pos; ) 

m 

rj void drawChartLine(Graphics g, String str_in, String str_type, double dbl_maxfreq) { 

H30 double dbl_x, dbl_inc, dbl_hmax; String str last, str next; 

str last str_next = ; dbl x = 25; gi count = 0; 

dbljimax = GetHMax(str_in); if(str type, equals IgnoreCase("S")) dbl hmax = 100; 

35 

dbl inc = ((((double)getSize().width - 50) / gi count) * ((dbl maxfreq - 200) / 800)); 
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for (int i cnt = 0; i_cnt < str_in.length(); ) ", i c{ 

nt-H-)if(str_in.regionMatches(i cnt, 0, 
2))if(str_last.length(> 0) { 
{ 

g.drawLine((int)dbl x, StringTolnt(str_Iast, dbl hmax), (int)(dbl x + dbl inc), 
StringTolnt(str_next, dbl hmax)); 
dbl x = dbl x + dbl inc; 

str last = strjiext; str next = ""; i cnt-w-; 

else { str next = str nextxoncat(str in.substring(i cnt, i cnt + 1)); 
) 

if(str_nextlength() > 0) 

g.drawLine((int)dbl x, StringTolnt(str_last, dbl hmax), (int)(dbl x + dbl-inc), 
StringTolnt(strnext, dbl_hmax)); ) 

public void paint(Graphics g) { double dbl x, dbl-y; 

g.setColor(Color.black); g.drawLine(0, 0, getSize().width, 0); g.drawLine(25, 
getSizeO-height - 25, getSize().width - 25, getSize().height - 25); g.drawLine(25, 25, 
25, getSize().height - 25); 

for(int i cnt = 0; i_cnt < 9; i cnt++) { 

dbl_x = 25 + ((double)i cnt * (((double)getSize(). width - 50) / 8)); 
g.drawLine((int)dbl x, getSize(). height - 25, (int)dbl x, getSize(). height - 20); 
g.drawString(String.valueOf((i cnt * 100) + 200), (int)dbl x - 8, getSize(). height - 5); ) 

for(int i cnt = 0; i_cnt < 1 1 ; i_cnt++) { 

dbl-y = 25 + ((double)i cnt * (((double)getSize().height - 50) / 10)); 
g.drawLine(20, (int)dbl-y, 25, (int)dbl-y); 

g.drawString(String.valueOf(100 - (i cnt * 10)), 1, (int)dbl-y + 5); ) 
g.drawString("Data Overlay", (getSize().width / 2) - 30, 12); 

g.setColor(Color.red); drawChartLine(g, getParameter("Data 1 "), getParameter("Type 
1 "), Double. valueOf(getParameter("Max Freq 1 ")). double Value()); 
g.setColor(Color.blue); drawChartLine(g, getParameter("Data 2"), 
getParameter("Type 2"), Double.valueOf(getParameter("Max Freq 
2")).doubleValue()); g.setColor(Color. green); drawChartLine(g, getParameter("Data 
3"), getParameter("Type 3"), Double.valueOf(getParameter("Max Freq 
3")).doubleValue()); ) ) 
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